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Abstract 

We demonstrate the method to calculate the spatial distributions of the spin current and accu- 
mulation in the multi-terminal ferromagnetic/nonmagnetic hybrid structure using an approximate 
electro-transmission line. The analyses based on the obtained equation yield the results in good 
agreement with the experimental ones. This implies that the method allows us to determine the 
spin diffusion length of additionally connected electrically floating wire from the reduction of the 
spin signal. 



I. INTRODUCTION 



Spin-dependent electron transport phenomena in nanostructured ferromagnet(F) /normal 
metal(N) hybrid systems show intriguing characteristics in association with the spin transfer 
and accumulation. The devices based on spin transport phenomena have great advantages 
over the conventional electronic devices because of additional spin functionalities.— There- 
fore, understanding responsible physics on the spin diffusion process becomes essential to 
realize spin-electronic devices. Especially, electrical spin injection from F into N, semicon- 
ductor or superconductor is an important technique for developing such devices.— When a 
spin-polarized current flows across the junction between F and N, the spin-splitting in the 
chemical potential is induced in the N layer because of the sudden change in spin depen- 
dent electrical conductivity.^ This leads to the spin accumulation in the vicinity of the F/N 
interface. The injected spin can be flipped due to spin-orbit interaction, magnon, phonon 
scattering and etc. The length scale over which the traveling electron spin memorizes the 
initial direction is known as the spin diffusion length, an important measure to realize the 
efficient spin injection. 

First electrical spin injection into N was demonstrated by Johnson and Silsbee using a 
single-crystalline Al bulk bar.— They found that the Al bar has a long spin diffusion length 
of a few hundred microns. However, the obtained output spin signal was quite small in 
the range of nano ohms due to large sample dimensions. Current perpendicular to plane 
(CPP) giant magnetoresistance measurements in magnetic multilayers of alternating F and 
N layers enable us to perform more detailed study on the spin-diffusion process in the N 
layer because the relevant scale of the spin-diffusion length can be controlled by the spacer 
thickness.— This technique is suitable for the N layer with magnetic or even non-magnetic 
impurities which reduce the spin diffusion length dramatically— However, for Ns in which 
long spin diffusion lengths are expected, the precise estimation of the spin diffusion length is 
unattainable because of difficulty in preparing CPP device with the thickness of the N spacer 
as thick as the spin diffusion length. On the contrary, planar mesoscopic spin-valve device 
is suitable to study the spin dependent electron transport of the material which has long 
spin diffusion length. Recently, Jedema et al. succeeded in detecting the clear spin signal 
in the lateral structure even at room temperature by using nonlocal spin-valve (NLSV) 
measurements similar to the Johonson's potentio-metric method.— In order to detect clear 
spin signals in the NLSV measurements, the spacing between the F injector and detector 
should be shorter than the spin-diffusion length of the material. This means that, for the Ns 



with spin diffusion length of sub hundred nanometers, the NLSV technique cannot be applied 
because of the technical limit of device fabrication. We have experimentally demonstrated 
that in such lateral structures an additional floating probe which does not carry the charge 
current, affects the spin transport when the probe is located within the spin diffusion length 
from the spin injector.— In this article, we further extend the experiment for quantitative 
estimation of the spin-diffusion length using the spin absorption effect, (i.e. spin sink) 

II. BASIC EQUATIONS 

The continuous charge current Ic in the F subsection provides the spin current 1$ = en pic, 
with cvf of the spin polarization in the F subsection. When the spin current flows across 
the F/N junction, the spin splitting in the chemical potential is induced in both F and N 
subsections, as shown in Fig. 1(a). Here, we devide the spin splitting chemical potential into 
two components. One is the voltage drop due to the charge current with the ohmic resistance 
in each layer shown in Fig. 1(b). In this case, the distribution of the chemical potential 
is obtained by considering the conventional electrical circuit consisting of a series jucntion 
with two different resistances. A constant spin current Is is induced by the continuous 
charge current in the F subsection. However, no spin current exists in the N subsection. 
The other one is the spin dependent the chemical potential without the contribution of 
ohmic-resistance voltage drop as shown in Fig. 1(c), where spin currents are generated both 
in the N and F layers. Since the sorce of the induced spin currents Is is originated from the 
charge current of the F layer, the relation Is = Isf + -^sn is satisfied at the interface. Here 
isF and -Zsn are respectively the spin currents in the F and N subsections induced by the 
spin splitting of the chemical potential. Thus, the F/N junction with the charge current Ic 
acts as a spin-current source with the magnitude of ctplc- 

We consider the spin-splitting component of the chemical potential in which there is no 
charge current. The essential features for the spin-current and spin- accumulation distri- 
butions in the diffusive regime are described by the spin-dependent Boltzman equation.— 
In general, from the one-dimensional diffusion equation, the induced spin-splitting voltage 
AVs = A/i/e = (/if — /ij,)/e, where /i| and /i; are the chemical potential of the up- and 
down- spins, respectively, is given by 

AVs = V + e~^ + V.e^. (1) 

Here A is the spin diffusion length. In the N layer, using the spin-dependent ohmic law 



= (— a^S/2)dV[ t i/dx, the spin current I$$ = — ^J. * s calculated as 

7 SN = _iLiL(^ +e ~^ - VCe^) = -— (V+e~^ - V_e^), (2) 
2An -Rsn 

where, Sn, <tn and An are the cross sectional area, the conductivity and the spin diffusion 
length of the N subsection, respectively. The spin resistance i?sN is a measure of the difficulty 
for spin mixing and is defined as 2An/(on<Sn)- The situation described by Eqs. 1 and 2 is 
equivalent to the electrical transmission line of the characteristic impedance i?sN with the 
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attenuation constant 1/An-- 

We can also obtain the similar relation for the F layer. In the F subsection, the spin 
current igp is given by 

a F S F A .dVt ,„ .dV,\ 
Isr = -y^ (a + a F )-^ - (1 - a F )^J (3) 

where, S F and ap are the cross sectional area and the conductivity in the F subsection, 
respectively. Since there is no charge current, the relation a F dAVs/dx = — {dV^/dx + 
dV±/dx) is satisfied. Therefore, the spin current I$ F in the F layer is given by 

Isf = (1 ~ ffi**^ (V + e~^ - K_e*) = -^(V + e-t - V.e^), (5) 

ZAp Ksp 

where, Ap is the spin diffusion length in the F subsection. Thus, in the F subsection, R SF 
is given by (1/(1 — a F ))(2X/(a F S F )) instead of 2An/(on<Sf)- ^qs. ^ and 5 mean that the 
characteristic spin resistance over the spin-diffusion length is comparable to the characteristic 
impedance for the spin current and spin-splitting voltage. 

Using these relations, we can simplify the calculation of the spatial distributions of spin 
current and accumulation in complex multi-terminal F/N hybrid structures consisting of 
different characteristic spin resistances. We will show that the characteristic spin resistance 
Rs is an important measure to design the spin-dependent-transport property of the system. 

First, we consider a cascaded transmission line of two different characteristic spin resis- 
tances as shown in Fig. 2(a). Here, the first and second subsections have the characteristic 
resistance Rsi with the spin-diffusion length Ai and the characteristic resistance Rs2 with 
the spin diffusion length A2, respectively. The length of the first subsection is d and the 
second subsection extends to infinity. The basic equations AVsi(a;) = Vi + e~ x l Xi + Vi_e x / Xl 
and Isi(x) = (Vi + e~ x ^ Xi — Vi~e x / Xi ) / Rsi are valid for each subsection. In the second line, the 
coefficient V2- of e x / X2 is zero because of its infinite length. Therefore, we obtain AV^ix) = 



y 2+e z/a 2 an( j / S2 (x) = AVs2(^)/-Rs2- This relation and the continuities of Is and AVs at 
the junction yield V 1+ = e d / x V 2+ (R S i + Rs2)/2R S 2 and Vi_ = e- d / A K 2+ ( J R S2 -R$i)/2R S2 . We 
should note that the situations where Vi_ > and Vi_ < respectively correspond to the 
spin-current absorption into the second line -and the reflection from it due to the resistance 
mismatch at the boundary. We then define the transmission coefficient T as the ratio of 
the spin splitting voltage Vjgi at the junction to the voltage Vso at the input terminal. The 
transmission coefficient T = AVjgi/AVso is given by 

T = ^2 = 9. ( 6 ) 

Rsi sinh(d/Ai) + R S 2 cosh(d/Ai) sinh(d/Ai) + Q cosh(d/Ai) ' 

where, Q is the ratio Rsz/Rsi- Defining the input spin resistance as AVso/iso where /so 
is the spin current at the input terminal, We obtain the input spin resistance in the series 
connection i?s_s er ies as 

^ _ i?sisinh((i/Ai) + i? S 2Cosh((i/Ai)^ _ sinh(d/Ai) + Q cosh(ci/Ai) ^ 

cnes Rsi cosh(d/ Ai) + i?s2 sinh(d/ Ai) cosh(c?/ Ai) + Q sinh(d/ Ai) 

The input spin resistance -Rs-ParaiM in a parallel connection shown in Fig. 2(b) is also 
calculated as 

„ RsiRs2 _ RS2 rs x 

-KS-Parallel - p D ~ \ . n - [°) 

ttSl + -KS2 J- + V 

Above basic equations can be applied for calculating the characteristics in multi-terminal 
spin-transport systems. 

In the above calculation, we assume a transparent interface between the different subsec- 
tions with zero resistance. This assumption is justified only in limited materials. However, 
in this article, we treat Cu/Py and Cu/Au interfaces. These interfaces are known to have 
small resistance compared with other interfaces such as Cu/Pd and Cu/Pt .^ 11 Therefore, 
we neglected the interface resistance. 



III. RESULTS AND DISCUSSIONS 



A. Single junction 

We now consider a single F/N cross junction with the charge current Ic flowing across 
the F/N interface as shown in Fig. 2(c). As mentioned above, the charge current across the 
junction plays a role of the spin current source with the magnitude of a.plc- The induced 
spin current diffuses into each subsection according to the magnitude of characteristic spin 
resistance. This can be described as the parallel spin resistance circuit consisting of two F 



and two N subsections shown in Fig. 2(c). We can now calculate the induced spin-splitting 
voltage Vso at the interface as 

AV S0 = a F I c , y (9) 

A-^SF + -KSNj 

The spin splitting decays exponentially over the spin diffusion length in each subsection. 



B. Double junction 

Unfortunately, we cannot experimentally determine the spin splitting of the chemical 
potential in the above single junction although the spin splitting can be induced in the N 
subsection. When an additional F probe (F2) is connected to the N subsection within the spin 
diffusion length from the spin-polarized current injector as shown in Fig. 2(d), the voltage 
difference between the F 2 and N subsections can be detected as the boundary resistance 
caused by the spin splitting. Remarked here is that the spin-current and chemical-potential 
distributions are significantly affected by the additional wire. We discuss the double F/N 
junctions on the basis of the transmission model. In this case, the input spin resistance 
Rsi for the nonlocal lead consisting of the N and F layers changes from -Rsn, resulting in 
the different spin-splitting voltage from the single junction. The Rsi can be calculated by 
considering a cascaded transmission line of i?sN and i?sd with the junction at the position 
d. Here, i?sd is the spin resistance from the detector junction and is given by 

-Rsn-Rsf/2 -Rsf /ir .N 

^ Sd " w^sn - 2Tg- (10) 

From Eq. i?si is deduced as 

i? S N smh(d/\ N ) + R S dCosh.(d/\ N ) 2i? S N ni >. 

R SN cosh(d/ A N ) + -Rsd smh(d/ A N ) 1 + e 1 ^ (1 + 2Q) 

The total spin resistance of the device Rs is given as the parallel sum of i?N, two Rf and 
Rsi '■ 

l 2 1 Y 1 QRsN\e^Q + 2smh(d/X N ^j 



R * = \ir + ir + ir) = ~n^ — (12) 

SN ^sf Rsi/ 2(e^Q(2 + g) + 2sinh(rf/A r 



Also using Eq. |Hl the transmission coefficient of the device T2 corresponding to the ratio of 
the spin-splitting voltage between the injector and detector F/N interfaces can be calculated 

as 

rp 2 = ^Sd = Q , 13 , 

2 i? S N sinh(d/ A N ) + R Sd cosh(rf/ A N ) e ^Q + 2 sinh(d/A N ) 



The induced spin-splitting voltage at the injector junction is given by IsRs- Then, the 
induced spin-splitting voltage AVs2 at the detector junction can be calculated by AVs2 = 
T 2 IsRs- Since the NLSV signal AV N2 is given by opAVsir^) we obtain 

AVn2 = x aFhQ2R ^ . (14) 

2e x N (2 + Q) +4sinh(rf/A7v) 

We experimentally demonstrate the precise estimation of the spin diffusion lengths of F 
and N strips and the spin polarization of F strips using Eq. El The lateral spin valves for 
this study consist of double Py/Cu junctions. Figure 3(a) shows a SEM image of the typical 
fabricated device. Here, the width of both Py and Cu strips is 100 nm and the thicknesses 
of Py and Cu strips are 30 nm and 80 nm, respectively. The resistivities of the Py and 
Cu strips are respectively 26.8 fiQcm and 2.08 fiQcm at room temperature. We vary the 
distance between the injector and detector junctions from 270 nm to 700 nm and measure 
the NLSV signal at room temperature as a function of the distance. In order to compare 
quantitatively Eq. El with the experimental spin signal, the cross sections Scu and Sp y have 
to be known to obtain the spin resistance of each strip. We thus have to carefully estimate 
the cross sections. From previous reports, Ac u is of the order of several hundred nanometer 
and Ap y is a few nanometers . 7 ' 13 ' 14 ' 15 Since the spin-current flows along the Cu strip over a 
few hundred of nanometers the area Squ for estimating Rq u should be given by the cross 
section of the Cu strip (100 nm x 80nm). On the other hand, in the Py strip, the cross 
section is not appropriate for the area Sp y for -RsPy because the spin current diminishes in 
the vicinity of the junction. Therefore, for the Py strip Sp y should be the junction-area 100 
nm x 100 nm.— 

A NLSV signal at a distance of 270 nm as a function of the external magnetic field parallel 
to the Py strip is shown in the inset of Fig. 3(b). We observe the clear spin-accumulation 
signal at room temperature. Figure 3(b) shows the NLSV signal as a function of the distance 
d. The obtained signal decreases monotonically with increasing the distance d due to the 
spin relaxation. This result is fitted to Eq. El As shown in the Fig. 3(b), the fitted 
curve is in good agreement with the experimental results. From the fitting parameter, we 
obtain the spin-diffusion length of the Cu strip and that of the Py strip as 500 nm and 3 
nm, respectively. And also, the polarization «p y of Py strip is determined as 0.25 at room 
temperature. Here, Ac u = 500 nm is quite long compared with other reported values.- 7 ^ 4 - We 
notice that the resistivity of our Cu is rather smaller than that of the reported one. This 
supports the long spin diffusion length of our Cu strip. The Ap y = 3 nm and ap y = 0.25 are 
reasonable although the values are slightly smaller than other group's values. This may be 



due to rather higher resistivity of our Py than others. Using these values, we can calculate 
the characteristic spin resistance of each strip as -Rscu = 2.60 Q and -Rsp y = 0.34 fl. 



C. Triple junction 

We extend the above formalism to a triple junction in Fig. 2(e). In this case, the output 
spin-splitting voltage induced between the F2 and N strip by the current injection from 
Fl is influenced by the M strip which absorbs the induced spin current. The middle strip 
has the spin resistance of Rsm with the spin diffusion length Am- The magnitude of the 
spin-current and spin-accumulation can be calculated by using the method explained above 
. From the calculation, the position of the middle wire was found not to affects the spin 
signal. Therefore, for simplicity we assume that the middle wire is located in the center 
between Fl and F2 strips which have the same spin resistances of Rsf- When the current 
is injected from Fl, the NLSV voltage AVn 3 at junction of F2/N is given by 

2 (-2 (1 + Q) + (2 + Q (2 + Q) (1 + Q M )) cosh(^) + (2Q M + Q (2 + Q) (1 + Qu)) sinh(^ 

(15) 

where, Q and Qm are respectively defined as Rsf/Rstss and _Rsm/-Rsn- d is the distance 
between the injector and detector. Qm = oo corresponds the case without M strip. 

The lateral spin valves consisting of the triple junctions are then prepared with changing 
the material for the M strip. We know the spin resistances of Py and Cu strips from the 
previous double-junction experiment. Therefore, the Py and Cu wires were used as F and 
N strip, respectively. Cu, Py and Au are chosen as the material of the middle strip. The 
width and thickness of each strip are the same as those of the double junction. Figure 4 
shows a fabricated lateral spin valve of the triple junctions with a Au M strip. The distance 
d between injector and detector Py/Cu junctions is 650 nm. The M wire is located at an 
intermediate position between the injector and detector. The NLSV is measured with the 
probe configuration in Fig. 4(b). Using the spin resistances -Rscu = 2.60Q, Rsp y = 0.34Q 
and Eq. [THl we obtain the following equation between the detected spin signal -Rnlsv and 
-Rsm 

_ 7.02/WK]] 
Rsu[n] ~ 1.27-4.41iW[mfi]- ( 6) 

Figure 5(a) shows the NLSV signal of the device with the Cu M strip. We obtained the 

spin signal of 0.18 mfl smaller than the value of 0.27 mfl obtained without the middle wire. 

This is because a small amount of the spin current is absorbed into the middle Cu strip. 



Figure 5(b) shows the NLSV signal of the device with the Py M strip. Here the difference 
in the switching field between Figs. 5(a) and (b) due to the shape of both injector- and 
detector-F strips to the large pad.— The spin signal exhibits a drastic reduction to 0.05 mQ, 
much smaller than that without the M strip. This is also due to the spin-current absorption 
into the Py M strip. However, because of the small spin resistance of the Py strip, the spin 
current is preferably absorbed into the Py M strip, resulting in much smaller spin signal in 
the NLSV measurement. Figure 5(c) shows the NLSV signal of the device with the Au M 
strip. The obtained spin signal is 0.08 mfl and shows the large reduction of the spin signal 
similar to that of the M Py strip. This implies that the Au wire has smaller spin resistance. 

Figure 6 shows the spin resistance of the M strip R SM as a function of -Rnlsv based on Eq. 
ITol From the equation, we obtain the spin resistances of the Cu, Au and Py stirp as 2.67, 
0.62 and 0.33 Q. The experimental values for the Py and Cu M strips are quantitatively in 
good agreement with the ones obtained from the double junction experiment. This indicates 
that Eq. El is valid for estimating the spin resistance of the M strip. For the Au strip, since 
the resistivity of the Au strip is 5.24 /jficm, we obtain the spin-diffusion length of the Au 
strip as 60 nm. This is in good agreement with the other reports^ and proves that the Au 
strip has short spin-diffusion length due to the strong spin-orbit interaction.— Thus, using 
this method, we can estimate the characteristic spin resistance and spin diffusion length of 
the material. 



IV. CONCLUSION 



We have shown the calculation of the spin-current and spin-splitting- voltage distributions 
in the multi-terminal F/N hybrid structure by introducing the characteristic spin resistance 
and similarly treating the system as an electro-transmission line. The analyses based on this 
method give the values in good agreement with the experimental results. The spin diffusion 
length of our Cu strip was found to be 500 nm quite long at room temperature. The model 
is extended to the triple F/N junctions consisting of a conventional lateral spin valve with an 
additionally connected floating strip which is located in between the injector and detector. 
When the spin resistance of the additional strip is much smaller than that of the N strip, 
the spin currents tend to flow into the additional strip to be relaxed. We also emphasize 
that, using the absorption effect, the spin-diffusion length of the additional strip can be 
estimated quantitatively from the reduction of the spin signal in the NLSV measurement. 
The situation of the present experiment may be directly compared with the spin pumping 



effect observed in magnetic multilayers.— 
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FIG. 1: (a) Schematic illustrations of the spin splitting chemical potential induced by spin injection 
in F/N junciton. (b) A component of the voltage drop due to the ohmic resistance and (c) that of 
the spin splitting without charge current. 

FIG. 2: (a) Schematic illustrations of junction connected in series (a) and in parallel (b). with 
different characteristic spin resistances The distribution of the spin current and spin splitting can 
be similarly treated as an electrical transmission line. Schematic illustrations of (c) a single cross 
F/N junction with the equivalent circuit for the transmission calculation, (d) that of a double F/N 
junction and (e) that of triple junction. 

FIG. 3: (a) SEM images of the typical lateral spin device consisting of the double Py/Cu junctions, 
(b) NLSV signal measured at room temperature as a function of the distance d between F injector 
and F detector. The inset of (b) is a obtained NLSV curve for the distance d = 270 nm. 

FIG. 4: SEM images of the typical lateral spin device consisting of the double Py/Cu junctions 
and the middle Au/Cu junction. 

FIG. 5: (a) Non-local spin-valve curve with the Cu middle wire, (b) that with the Py middle wire 
and (c) with for the Au middle wire. 



FIG. 6: Calculated spin resistance as a function of the obtained non-local spin signal i?NLSV- 
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